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Abstract 

The article considers the successful implementation of relativistic equation-of-motion 
coupled cluster method for the electron attachment problem (EA-EOMCC) at the level 
of single- and double- excitation approximation. The implemented relativistic EA- 
EOMCC method is employed to calculate ionization potential values of alkali metal 
atoms (Li, Na, K, Rb, Cs, Er) and the vertical electron affinity values of LiX (X=H, E, 
Cl, Br), NaY (Y=H, E, Cl) starting from their closed-shell configuration. Both four- 
component and exact two component calculations are done for all the opted systems. 
Eurther, we have shown the effect of spin-orbit interaction considering the atomic sys¬ 
tems. The results of our atomic calculations are compared with the values from the 
NIST database and the results are found to be very accurate (< 1%). 
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1. Introduction 


A considerable growing interest is noticed in recent years in the study of negative 
ions as negative ions have significance in maiw areas of physics like in astrophysics, 

mm 


plasma physics and surface physics 


The electron affinity (EA) is an im¬ 


portant quantity of these ions. The precise measurement of EA of atomic or molecular 
systems is always a challenge as the resulting negative ion is difficult to handle. De¬ 
spite of the complexity in the measurement, there have been significant advances in 
the experimental techniques like laser photodetachment electron spectroscopy (LPES), 
laser photodetachment threshold spectroscopy (LETS), accelerator mass spectroscopy 
(AMS) and photodetachment microscopy, et cetera that are capable of precise measure¬ 


ments of EA of an atomic system 


oscOTwet ce 

iflflIaS]. 


However, the situation is inappreciative 


in achieving such an extent of accuracy in the molecular systems due to the possibility 
of structural change on attachment of an extra electron. Therefore, it is an outstanding 
challenge for the computational physicists to complement these atomic measurements 
as well as for new predictions for the future purpose. 

The computational prediction of EA is difficult due to the absence of long-range 
Coulomb field outside of a neutral precursor. Therefore, an extra electron is solely 


bound through correlation with other electrons 


iJQ. 


Moreover, most of the theo¬ 


retical calculations are based on the quantum chemical basis set methods. Thus, the 
finite size of the basis and unbalanced treatment of electron correlation in the atomic or 
molecular system and in the resulting ion are the sources of error | nl The attach¬ 
ment and detachment of an electron to a neutral species involves different forces. The 
attached extra electron to the neutral atom polarizes the electronic shell of the atom. 
As a result a dipolar electric field is generated which binds the extra electron with the 
other electrons. The charge distribution of the electron cloud, particularly the electron- 
electron correlation effects decides the stability of the negative ion. These interactions 
do not play much role in most of the neutral atoms as well as in positive ions where di¬ 
rect electrostatic force is the dominant factor for the stability of the neutral atom or the 
positive ion. On the other hand, these effects dominate in the negative ions. Therefore, 
the calculations of EA values of both atomic and molecular systems are challenging 
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and is a real test for the performance of a many-body method. It is an established fact 
that not only the electron correlation but also the effect of relativity play a definite role 
in accurate description of the eigenstates of heavy atomic and molecular systems ll^ . 
It is, therefore, in such a case a highly correlated many-body method, capable of simul¬ 
taneous treatment of relativity and electron correlation is required due to the intricate 
coupling between these two effects | 


QQQ]. 


The relativity has a greater role towards the core orbitals and practically important 
for all the elements. The effects of relativity are incorporated in the electronic structure 
calculations by the choice of the Hamiltonian. The consideration of Dirac-Coulomb- 
Breit Hamiltonian without the quantum electrodynamics effects (QED) is sufficient for 
most of the relativistic electronic structure calculations using four-component wave- 
function. However, in actual practice the Dirac-Coulomb Hamiltonain is most com¬ 
monly used where two-body Coulomb interaction operator is added to the Dirac Hamil¬ 
tonian {Ho)- Although, the form of the Coulomb operator is same as in the non- 
relativistic theory, however, the physical content is different as it takes care of the 
spin-same orbit interaction. This type of truncation in the two-body interaction does 
not effect much for most of the chemical purposes ll^ . However, for very accurate 
studies of molecular spectra including fine structure, the inclusion of spin-other-orbit 
interaction and spin-spin interaction are required which can be done with the full in¬ 
clusion of the Breit part of the two-body interaction. The relativistic calculations using 
four-component wave function are very expensive from the computational perspec¬ 
tive. A lot of effort has been made to simplify the equations. The calculation of the 
small component of the wave function is the most challenging part of the computation. 
If a basis set is expressed in terms of contracted Gaussian functions, then the num¬ 
ber of required primitive Gaussian functions for the small component is about twice 
the number of the large component with the imposition of the kinetic balance condi¬ 
tion. On the other hand, the small component has a very minor contribution in the 
calculated values; therefore, it makes sense to look for an approximation. There are 
a number of Hamiltonians in between the scalar non-relativistic and four-component 
relativistic ones. However, the inclusion of the spin-orbit interaction requires at least a 
two-component description, though it will essentially increase the computational cost 
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due to the appearance of complex algebra in place of real algebra. The electron cor¬ 
relation methods in the no-pair approximation require the transformation of the matri¬ 
ces from the atomic orbital (AO) basis to the molecular orbital (MO) basis. The spin 
coordinates of the electrons can be represented in terms of quaternion algebra in the 
four-index transformation step which helps to go from complex four-component to a 
two-component quaternion form. Therefore, the MO coefficients become quaternion 


and can be represented in terms of real matrices [[IS 




The generation of a two-component Hamiltonian from the parent four-component 
Hamiltonian is the most preferred choice for the purpose which includes the spin-orbit 
interaction with a lesser cost as compared to the four-component Hamiltonian. The cen¬ 
tral idea behind the generation of a two-component Hamiltonian is that it should repro¬ 
duce the positive-energy spectrum of the parent Hamiltonian. Foldy and Wouthuysen 
proposed an idea to decouple the large and small component by a unitary transforma¬ 
tion of the four-component Hamiltonian. Another well known approach is the elimi¬ 
nation of the small component from the wavefunction. However, these two approaches 


can be shown to be equivalent ll20l] . The exact two-component approach (X2C) in 


the two-component framework is one such approach to reduce the computational scal¬ 
ing which uses elimination of the small component from the parent four-component 
Hamiltonian. The detailed description of the X2 C appro ach including various other 


two-component methods can be found in Ref. 113 


21 


221 ]. 


Over the years, the equation-of-motion coupled-cluster (EOMCC) method gained 
popularity among correlation methods for the treatment of electron correlation due to 


its simplicity and elegance. The idea of EOMCC l23 


24 


25 


m 


27 


28 


29 


3C 


nil is 


conceptually very simple and it is operationally a two step process: (i) solution of cou¬ 
pled cluster problem with the N electron closed-shell determinant as reference and (ii) 
construction and diagonalization of the effective Hamiltonian matrix for the Eock-space 
sector of interest in the configuration space. It takes into account of both the dynamic 
and non dynamic part of the electron correlation. The exponential structure of the cou¬ 
pled cluster operator takes care of the dynamic part of the electron correlation and non 
dynamic part is included by means of diagonalization of the effective Hamiltonian ma¬ 
trix in the configurational space. The diagonalization of effective Hamiltonian, by and 
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large is associated with the multi-reference theories, whereas EOMCC works within a 
single reference description to tackle the complex multi-configurational wavefunction. 
Further, the relaxation effect, which has an important role in proper description of the 
eigenstates is also taken care. The multiple roots can be addressed in a single calcu¬ 
lation and each of the states are treated with equal weightage. The EOMCC method 
behaves properly at the non-interacting limit but not rigorously extensive (only for the 
core-core and core-valence interactions) due to the linear structure of the EOM operator 


1 3211 . The EOMCC is in close kinship with the coupled cluster linear response theory 
341] and symmetry adapted cluster expansion configuration interaction 


(CCLRT) I13J 


(SAC-CI) method 135 


361] . It is worth to note that the transition energy calculated 


using CCLRT is identical with the EOMCC method for the one valence problem but 
the transition moments is identical only when it is represented as a energy derivative 


in EOMCC framework. Chaudhuri et al l37 


3811 applied relativistic CCLRT for the 


ionization problem of atomic systems with spherical implementation. Besides these 
two methods, effective Hamiltonian variant of Lock space multi-reference (FSMRCC) 


theory l39 


40 


41 




43 


441] always comes in the discussion on EOMCC as these two 


methods produce identical results for the one valence problem. The amplitudes of all 
the lower sector including the sector of interest are involved in the FSMRCC theory. On 
the other hand, EOMCC deals with the amplitudes of the (0,0) sector and the sector of 
interest. Therefore, both the approaches are eventually produce the same result for the 
one electron attachment or detachment problem. The EOMCC is free from the problem 
of intruder due to its Cl (configuration interaction method) like structure, which is as¬ 
sociated with the effective Hamiltonian variant of the FSMRCC theory. There are ways 
in the FSMRCC framework to handle the problem of intruder such as the eigenvalue 
independent partitioning technique of Mukherjee (EIP-FSMRCC) 


321 145[] and the in¬ 


termediate Hamiltonian variant of the FSMRCC (IH-FSMRCC) theory , 


47 


4 ^]. 


Recently, Blundell implemented relativistic EOMCC method for the electron affin¬ 
ity problem and applied to calculate fine-structure splittings in high-lying states of 
rubidium atom 1491] . The implemented version of Blundell is applicable only for the 
purpose of atomic calculations as they have used the spherical implementation which 
allows the separation of radial and angular parts. Therefore, the evaluation of radial in- 
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tegrals is only required and the angular part will add up to it as a multiplier. The radial 
integrals can be evaluated numerically. Such a separation is not possible in molecular 
systems due to the absence of spherical symmetry. In our implemented version, we 
have used the one- and two- body matrix elements, which are evaluated in the Carte¬ 
sian coordinate system. The Cartesian coordinate system does not allow one to exploit 
the spherical symmetry to separate the matrix elements into radial and angular parts. 
Furthermore, the anti-symmetrized two-body matrix elements are used in this coor¬ 
dinate system calculations, which is not possible in the spherical implementations as 
angular factor will be different for the direct and exchange part of the two body matrix 
element. Thus, our implemented version is a general one, applicable to both atoms 
as well as molecules starting from their closed-shell reference state configuration. It 
should be noted that the spherical implementation is much more complex than that of 
the molecular calculations, but it is favorable from the computational point of view as 
it requires only the solution of radial integrals. Therefore, atomic calculations are com¬ 
putationally easy, which allows to correlate more number of electrons and amenable to 
use huge basis for the correlation calculation to achieve a better accuracy. 

In our recent work, the performance of the fully four-component EOMCC has al¬ 
ready been established for both atomic and molecular systems for the single ioniza- 


|3QQ- 


Therefore, in this article, we focus 


tion and double ionization problem 
on the implementation of relativistic EOMCC method for the electron affinity prob¬ 
lem applicable to both atomic and molecular systems. The implemented EA-EOMCC 
method is employed to calculate ionization potential of open-shell atomic systems start¬ 
ing from their singly positive closed-shell configuration. Eurther, the vertical EA values 
of molecular systems are also calculated. Both four-component and exact two compo¬ 
nent (X2C) calculations are done for all the considered systems. The effect of spin-orbit 
interaction is shown for the atomic systems in the EOMCC framework. 

The manuscript is organized as follows. The EOMCC theory in regard to the elec¬ 
tron attachment problem is briefly described in Sec. |2]and the computational details of 
our calculations are presented in Sec. [3 We have discussed our results in Sec. |4]and 
finally made concluding remarks in Sec. |5] We are consistent with the atomic unit if 
not stated explicitly. 
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2. Theory 


In the EOMCC method the target excited state of single electron attached state 
is defined as 

|vI/fc) = <+i|vI/o),fc = l,2,... (1) 

Here, the is a linear operator, which on acting on the single reference coupled 

cluster (SRCC) ground state wave function |4'o), generates the excited state wave 
function The operator takes the form in the coupled cluster single-double 

(CCSD) approximation as 

+1 = i?i -P i?2 

= ^ ^ ^ ( 2 ) 

a b<a j 

The i?i and i ?2 operator are diagrammatically represented in Fig. [T] The Ri is a one 
particle (Ip) creation operator and R 2 is a two-particle and one-hole (2p-lh) creation 
operator. The circled arrow is just to represent that overall it is a one electron attach¬ 
ment process. 



Figure 1: Diagrammatic representation of Ri and R 2 operator. 


The Schrodinger equation for the ground state (k=0) is 

= AFiol^-o), (3) 

The electron attached states (k=l,2,...) is written as 

HNRkl'i’o) = AEkRkl'i’o) (4) 

The above equation on multiplication with a non-singular operator (where T is 
the coupled cluster excitation operator) in the course it is assumed that Rk commute 
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with T (as strings of same quasi-particle creation operator) with some mathematical 
manipulation leads to equation of motion with respect to the Rk operator, 

vfc. (5) 


In the above equation, /S.Ek is the energy change associated with the electron attach¬ 
ment process and Hn = He^ — (<l>o|e~^'f^e^|‘&o) is the similarity transformed 
normal ordered effective Hamiltonian. In our case it is the Dirac-Coulomb Hamilto¬ 
nian, which is given by 



where and /3 are the usual Dirac matrices, ViA is the nuclear potential and — is the 
electron-electron repulsion potential. The orbital energies are scaled with respect to the 
free electron rest mass energy {(?), which is zero in the non-relativistic case. We have 
chosen a correlated determinantal space of |<I>“) and |<I>“^) (Ip and 2p-lh) with respect 
to the Dirac-Hartree-Fock determinant (|$o)) to project the above equation to get the 
desired electron affinity values, 


($“[iJAr,Efc]|$o) = AFifci?^, (7) 

($f [Fat, i?fc]|$o) = (8) 


In Figs. |2]and[3 the contributing diagrams for the Ip and 2p-Ih space is presented, 
respectively. The evaluation of these diagrams requires the solution of the coupled- 
cluster amplitude equations. The converged Ti and T 2 amplitudes are contracted with 
the two-body matrix elements of the Hamiltonian matrix to construct one body, two 
body and three body intermediate diagrams. The intermediate diagrams are categorized 
into fpp, fhp, fhh, Vhppp, Vpppp, Vphph, Vppph and W. Here /’s, t/’s and W stands for 
one-body, two-body, and three-body intermediates, respectively. We have followed a 


recursive intermediate factorization scheme as described in Ref. 11531] to evaluate these 


intermediate diagrams. The factorization scheme in the construction of intermediate 
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diagrams saves enormous computational resources. The matrix elements correspond¬ 
ing to the three-body intermediate diagram are not stored rather computed on the fly. 
The programmable algebraic expression for the diagrams corresponding to the projec¬ 
tion of Hamiltonian to Ip and 2p-lh excited determinantal space are presented in Eq. |9] 
and[T0l respectively. We have used the standard notation (/(out, in)) and {V (left out, 
right out, left in, right in)) for one-body and two-body intermediate matrix element. 
In Eqnl9]and[T0l i,j,... {a,b...) stands for hole (particle) index. P is a permutation 
operator and any odd permutation introduces a negative sign. The Davidson algorithm 
1 5411 is implemented to get the desired set of eigenvalues and eigenvectors. 


APfcP“ = ^ /pp(a, b)r^ -f ^ fhp{j, byf + 0.5 ^ Vhppp{j, a, b, c)rf V a (9) 

b j^b j,b^c 


APfcPf = Piab) fPPib, c)r“ - ^ hh{j, + 0.5 ^ Epppp(a, b, c, d)rf 

c j c,d 

-P{ab) E ^phphi^i ^ ^ ^ppph (Q; J 

c,k c 

-0.5 ^ Vhhppik,l,c,d)rlX^ V (/ b<a) (10) 

k,l,c,d 
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a 

(> 




Figure 2: Diagrams contributing to the Ip block. 




Figure 3: Diagrams contributing to the 2p-lh block. 
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Table 1: SCF and correlation energy from MBPT(2) and CCSD of alkali metal ions. 




X2C 



4C 


Atom 

SCF 

MBPT(2) 

CCSD 

SCF 

MBPT(2) 

CCSD 

Li+ 

-7.237045 

-0.038520 

-0.042313 

-7.237174 

-0.038490 

-0.042284 

Na+ 

-161.885871 

-0.352072 

-0.355824 

-161.895637 

-0.351504 

-0.355258 

K+ 

-601.317915 

-0.655371 

-0.669654 

-601.378197 

-0.654132 

-0.668419 

Rb+ 

-2979.125369 

-1.603877 

-1.544953 

-2979.693217 

-1.600533 

-1.541667 

Cs+ 

-7784.579785 

-1.746204 

-1.654556 

-7786.643511 

-1.740037 

-1.648569 

Fr+ 

-24296.910671 

-1.635291 

-1.468312 

-24308.061505 

-1.634768 

-1.467783 








Table 2: Bond length(in A), SCF {E^^} and correlation energy from MBPT(2) and CCSD of LiX(X=H, F, Cl, Br) and NaY(H, F, Cl). 


Molecule 



X2C 



4C 


Bond Length [551 

SCF 

MBPT(2) 

CCSD 

SCF 

MBPT(2) 

CCSD 

LiH 

1.5957 

-7.987794 

-0.069223 

-0.080955 

-7.987928 

-0.069191 

-0.080923 

LiF 

1.5939 

-107.078887 

-0.400341 

-0.399171 

-107.084024 

-0.400097 

-0.398927 

LiCl 

2.0207 

-468.469738 

-0.606016 

-0.622235 

-468.511539 

-0.605515 

-0.621738 

LiBr 

2.1704 

-2612.134924 

-1.465099 

-1.404927 

-2612.603993 

-1.463721 

-1.403588 

NaH 

1.8874 

-162.602176 

-0.385821 

-0.397925 

-162.611952 

-0.385251 

-0.397357 

NaF 

1.9259 

-261.676902 

-0.694456 

-0.691459 

-261.691677 

-0.694011 

-0.691018 

NaCl 

2.3608 

-623.077913 

-0.846128 

-0.861815 

-623.129354 

-0.845488 

-0.861179 








3. Computational Details 


The DIRAC program package 156 , 


5711 is used to evaluate the required one-body 


and two-body matrix elements for the correlation calculation. Both the X2C and four- 
component calculations are done by using uncontracted finite atomic basis, which is 
consists of scalar real gaussian functions. The small component of the basis is linked 
with the large component of the basis through the restricted kinetic balance (RKB) 
condition. The RKB condition represents the kinetic energy properly in the non rela¬ 


tivistic limit and avoids the variational collapse 11581] . This is achieved by pre-projecting 


in scalar basis and unphysical solutions are removed by diagonalizing the free particle 
Hamiltonian. The RKB condition generates the positronic solution and electronic so¬ 
lution in 1:1 manner. The DIRAC program package uses Gaussian distribution nuclear 
model to take care of the finite size of the nucleus. The used nuclear parameters are 


taken as default values from DIRAC package 11591] . We adopted aug-cc-pCVQZ basis 


for Li+ ll60l] and Na+ ]l6ll] atom and all the generated orbitals are taken into consider¬ 
ation for the correlation calculations. Dyall.cv4z 11621] basis is opted for K+ and Rb+. 
We have neglected the virtual orbitals those energy is more than 500 a.u. for the K+ 
and Rb+ atom. The Cs+ and Fr+ are calculated using dyall.cvSz basis 11621] . The cut¬ 
off used for Cs+ atom is 1000 a.u. whereas for Fr+ atom, we have taken the orbitals 
having energy in between —25 a.u. to 100 a.u. in our correlation calculations. In the 
molecular calculations of LiF, LiCl, LiBr, we have chosen aug-cc-pCVTZ basis for Li 


atom ll60ll and cutoff of 100 a.u. for the virtual orbitals. The calculations of F and Cl 


are done using aug-cc-pCVQZ ll60ll6lll basis and for Br, dyall.cv4z 116411 basis is used. 


In LiH we have chosen aug-cc-pCVQZ basis 06011 for Li and aug-cc-pVTZ ll60ll for the 
H atom and none of the electrons are frozen for the correlation calculations. Aug-cc- 


pCVTZ basis is opted for both Na mill and Cl 06511 in the calculations of NaCl and a 
cutoff of 100 a.u. is used for the virtual orbitals. The single particle orbitals and two- 
body matrix elements are generated by taking account of C 21 , symmetry. Both X2C 
and four-component calculations of Rb and LiBr are done with the DIRAC 14 version 
and rest of the calculations are done using DIRAC 10. The implemented version of 
X2C SCF B in DIRAC 10 is capable of taking up to g harmonics but the opted basis 
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Table 3: Convergence pattern of electron affinity (in eV) of the C 2 (R= 1.243 A, Ref 15511 1 as a function of 
basis set. 


Basis 

Spinor 

Electron affinity 

Expt.[^ 

cc-pVDZ 

88 

2.6494 


cc-pCVDZ 

104 

2.6645 


aug-cc-pVDZ 

124 

3.1896 


aug-cc-pCVDZ 

140 

3.1949 


cc-pVTZ 

152 

3.1140 


cc-pCVTZ 

192 

3.1285 


aug-cc-pVTZ 

216 

3.3316 

3.30±0.1 

cc-pVQZ 

252 

3.2743 


aug-cc-pCVTZ 

256 

3.3412 


cc-pCVQZ 

316 

3.2851 


aug-cc-pCVQZ 

416 

3.3840 


aug-cc-pCVQZ“ 

488 

3.3853 



“ All the virtual orbitals are used for the EOMCC calculation. 


for Rb and LiBr require up to h harmonics to express the large component of the wave 
function. Therefore, these two calculations are done using DIRAC14 version. We have 
fixed a cutoff of 10“^^ to store the matrix elements for the intermediate diagrams as 
two-body matrix elements contributed negligibly beyond this limit. The convergence 
of 10“® is fixed for the solution of SRCC amplitude equations and 10“® for the David¬ 
son algorithm in the EOMCC part. A direct inversion in the iterative subspace (DIIS) 
of 6 is used in the solution of ground state amplitudes for all the calculations. The 
newly implemented relativistic EOMCC code is tested by comparing EA-EOMCC re¬ 
sults with the (1,0) sector ESMRCC code implemented in the DIRAC package as these 
two theories are supposed to produce identical results for one electron attachment pro¬ 
cess. The MBPT(2) correlation energy is identical whereas CCSD correlation energy 
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Table 4: Ionization potential values (in eV) of alkali metal atoms. 


Atom 

X2C 

4C 

NIST 1671 

Li 

5.3894 

5.3895 

5.3917 

Na 

5.1104 

5.1106 

5.1391 

K 

4.3419 

4.3423 

4.3407 

Rb 

4.1750 

4.1756 

4.1771 

Cs 

3.8861 

3.8872 

3.8939 

Er 

4.0579 

4.0603 

4.0727 


and the EA value are matching upto ten-digit and eight-digit, respectively. The test is 
performed with identical convergence cut off, equal number of DIIS space and without 
any cutoff in the intermediate diagrams. We have done the test over a series of atoms 
and molecules with various basis sets and successful in achieving similar agreement 
for all the considered system, independent of the choice of basis set. 

4. Results and discussion 

We have reserved this section of our manuscript to present numerical results of our 
calculations and to interpret the outcome of these calculations. The four-component 
and exact two component (X2C) EOMCC calculations are performed for all the con¬ 
sidered atomic and molecular systems starting from their closed-shell configuration. 

In Tables [Hand in we present numerical results of our SCE and correlation energy 
from MBPT(2) and CCSD calculation of singly positive alkali metal atomic systems 
(Li+, Na+, K+, Rb+, Cs+, Er+) and molecular systems (LiX (X=H, E, Cl, Br) and NaY 
(Y=H, E, Cl)) in their closed-shell configuration. The bond length of the molecular 
systems is also compiled in Table [H Eor atomic systems (Table [TJ, we noticed that 
the difference between MBPT(2) and CCSD correlation energies for both X2C and 
four component calculation keep on increasing as we go down the group. This trend is 
expected as the effect of correlation increases as the number of electron in the system 
increases. 
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Table 5: Relaive energy difference (in eV) of energy levels of atoms 


Atom 

2 S 


'Pl/2 



'P 3/2 



"P 3/2 -" 

Pl/2 

X2C 

4C 

NIST 

X2C 

4C 

NIST [671 

X2C 

4C 

NIST 

Li 

0.0000 

1.8494 

1.8495 

1.8478 

1.8495 

1.8495 

1.8479 

0.0001 

0.0000 

0.0001 

Na 

0.0000 

2.0837 

2.0840 

2.1023 

2.0859 

2.0862 

2.1044 

0.0022 

0.0022 

0.0021 

K 

0.0000 

1.6128 

1.6132 

1.6100 

1.6203 

1.6206 

1.6171 

0.0075 

0.0074 

0.0071 

Rb 

0.0000 

1.5610 

1.5616 

1.5596 

1.5909 

1.5915 

1.5890 

0.0299 

0.0299 

0.0294 

Cs 

0.0000 

1.3912 

1.3921 

1.3859 

1.4600 

1.4609 

1.4546 

0.0688 

0.0688 

0.0687 

Fr 

0.0000 

1.5198 

1.5220 

1.5172 

1.7276 

1.7298 

1.7264 

0.2078 

0.2078 

0.2092 








We have done a series of calculations to understand how the electron affinity value 
changes with the nature of basis set and cutoff in the orbital energies can be used 
without losing considerable amount of accuracy as EOMCC calculations for the EA 
problem are computationally costly. We have chosen C 2 as an example system for 
which experimental vertical EA value is reported in the literature. We have started our 
calculation with cc-pVDZ, which is a very small basis as it generates only 88 spinor 
for the beyond SCE calculations using a cutoff of 100 a.u. for the virtual orbitals. A 
few more calculations are also done by improving the nature of the basis functions. 
The EA value as well as the number of generated spinor in different basis are tabulated 
in Table [3 We have taken 1.243 A as bond length for the C 2 molecule, which is the 


experimentally reported bond length 0551] . The reported experimental value is 3.30±0.1 


eV 11631] . whereas our calculation yields 3.3840 eV for aug-cc-pCVQZ basis with a 


cutoff of 100 a.u. in the virtual orbital energy. On the other hand, without using 
any cut off, the result is 3.3853 eV. Therefore, a cutoff of 100 a.u. for virtual orbital 
energies and similar basis set can be used without loosing much accuracy for all other 
calculations to achieve a good agreement with the experiment. It will save enormous 
computational time without losing a significant amount of accuracy as contribution 
from the high energy virtual orbitals is very less in the correlation calculations. The 
reported experimental uncertainty is in the first digit after the decimal point. Therefore, 
it is hard to comment on the accuracy of our calculated results. It can be said that 
our results are also spanning same range starting from a reasonable basis. In Table 
m we report the calculated ionization potential values of atomic systems using both 
X2C and four-component EA-EOMCC method. We have started our calculations from 
singly positive alkali metal ions and applied EA-EOMCC method. The negative of the 
computed values are reported as ionization potential values of the open-shell atomic 
systems. These computed ionization potential values are compared with the values 
from the National Institute of Standards and Technology (NIST) database. A nice 
agreement with NIST values is achieved for all the considered systems. The maximum 
deviation is obtained for the Na atom, which is also in the accuracy of ~ 0.6% with 
the NIST value. The difference between the X2C with four-component results is in the 
fifth digit after the decimal for Li atom whereas the difference is about 0.01 eV for Er 
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Table 6: Vertical EA (in eV) values of LiX (X=H, F, Cl, Br) and NaY (Y=H, F, Cl). 


Molecule 

X2C 

4C 

Others 1681 

LiH 

0.2968 

0.2968 

0.247 

LiE 

0.3550 

0.3550 

0.340 

LiCl 

0.5526 

0.5526 

0.551 

LiBr 

0.6148 

0.6148 


NaH 

0.3218 

0.3217 

0.319 

NaE 

0.4848 

0.4849 

0.485 

NaCl 

0.6726 

0.6727 

0.672 


atom. 


We have calculated electron attachment energy to the pi /2 P 3/2 orbitals of the 
atomic systems in the X2C and four-component EOMCC framework and presented as 
a relative energy difference with respect to the '^S states. These results are compared 
with the values from NIST database and presented in the Table |5] The results of our 
calculations are found to be very accurate with respect to the NIST values. It is interest¬ 
ing to note that the energy gap between the pi /2 and p ^/2 is negligible for the Li atom 
but it keeps on increasing as the system size become havier. This can be explained by 
the fact that the effect of relativity increases with the increase in atomic number and 
thus, the spin-orbit coupling plays a significant role in heavier systems. 

In Table | 6 l we present the results of our calculations of vertical EA of molecular 
systems using both X2C EA-EOMCC and four-component EA-EOMCC method. Eur- 
ther, we have compared our result with the theoretical results calculated by Gutsev et 


al 11681] . They also employed EA-EOMCC method for correlation treatment to calcu¬ 


late the vertical EA values of the molecular systems. In their calculation, Gutsev et al 
misses the effect of relativity, which is included in our calculation in its four-component 
formalism. We have achieved a nice agreement for the atomic results and also for the 
vertical EA value of C 2 molecule. Therefore, it can be said that our calculated results 
for the molecular systems are also quite accurate though there is no reliable experimen- 
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tal data or any other values calculated using any variant of relativistic coupled cluster 
theory to compare with. However, the accuracy of the molecular calculations will not 
be that much accurate as compared to the atomic results. The reason behind this is due 
to the possibility of structural change on attachment of an extra electron to the neutral 
molecule depending on the polarity of the molecule. 

5. Conclusion 

The relativistic EOMCC method for the electron attachment problem applicable to 
both atomic and molecular systems is successfully implemented. To test the perfor¬ 
mance of the EA-EOMCC method, we applied to calculate ionization potential values 
of alkali metal atoms starting from closed-shell configuration. We have compared our 
calculated ionization potential values with the values from NIST database. We are 
successful in achieving less than 1% agreement with the NIST values. We have also 
presented molecular EA values of LiX (X=H, E, Cl, Br) and NaY (Y=H, E, Cl) using 
our relativistic EOMCC methods. 
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